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ABSTRACT 

The electromagnetic emissivity from QCD media away from equihbrium is studied in the framework of closed 
time path thermal field theory. For the dilepton rate a nonequilibriuni mesonic medium is considered applying 
finite temperature perturbation theory for ir — p interactions. The dilepton rate is derived up to the order g^. 
For hard photon production a quark gluon plasma is assumed and calculations are performed in leading order 
in the strong coupling constant. These examples are chosen in order to investigate the role of possible pinch 
terms in boson and in fermion propagators, respectively. The implications of the results for phenomenology are 
also discussed. 



' Contribution to the XXXVII-th Cracow School of Theoretical Physics, Dynamics of Strong Interactions, to appear in 
Acta.Phys.Pol. 



I. INTRODUCTION 



Dilepton and photon production is one of the interesting tools to study coUective effects in strongly interacting 
matter produced in ultrarelativistic heavy ion collisions . This is particularly evident from the recent CERN 
experimental data at SPS collision energy measured by the collaborations CERES/NA45 |] and HELIOS/3 
j^, which show qualitatively different behaviour of the dilepton yields in A- A as compared to p-p collisions. 
Extensive theoretical studies done over the last years ]7|-p^ suggest that this new behaviour of the measured 
dilepton yield could be related with the thermalization of the medium created in heavy ion collisions. Also the 
properties of the measured particle spectra at AGS and SPS energy seem to be well explained by the production 
from the thermal source. 

Thermalization of a hadronic medium should be a fast process. Recent calculations in kinetic theory show that 
only a few elastic particle collisions are already sufhcient to maintain thermal particle spectra, The chemical 
equilibration, in contrast, could be much slower as it requires a detailed balance between different reactions 
with particle number changing processes. Also recent predictions of QCD inspired models like HIJIING or 
Self-Screened Parton Cascade jl^ lead to similar conclusions for the partonic medium. These models indicate 
that the partonic medium created initially in ultrarelativistic heavy ion collisions appeared after very short 
time in local thermal equilibrium, however, chemical saturation is most likely never achieved before the system 
hadronises. 

The above discussions thus show that when analysing the properties of strongly interacting matter as produced 
in heavy ion collisions we have necessarily to take into account nonequilibrium effects. 

In this work we show how nonequilibrium effects could possibly influence the electromagnetic emissivity from 
QCD media. In particular we discuss soft dilepton production from tt — p interactions in an off-equilibrium 
mesonic medium and hard thermal photon production from an undersaturated quark gluon plasma. 

We apply the real time formalism of field theory with particular emphasis on the role and structure of the 
pinch singularities [|T^,Q which appear in propagators of particles in a nonequilibrium medium. 

The paper is organized as follows: in Section. 2 we introduce the Keldysh representation and discuss the 
structure of particle propagators in nonequilibrium media. In Section. 3 the dilepton rate from pi-rho scattering 
is derived and discussed. In Section. 4 we show that there is self-consistent dynamical screening which leads to 
a finite photon rate in an off-equilibrium quark gluon plasma. 



II. REAL TIME FORMALISM AND PARTICLE PROPAGATORS 



In order to see how nonequilibrium effects modify the structure of the particle propagators we first summarize 
some the important relations of the real time formulation of finite temperature field theory. 

In the context of the closed time path or the Keldysh representation the bare propagator in equilibrium at 
finite T acquires a 2 x 2 matrix structure of the following general form: [p^|-pl| 

where k = {kg, k), k"^ ~ k'^ — k"^ \ 9 denotes the step function, m is the T = mass and ri(fco) is the thermal 
distribution function 

"(fco) = . i ^ , , (11.2) 



where the upper (lower) sign in the denominator of the above equation (similarly as in (II. 1)) corresponds to 
bosonic (fermionic ) fields. The Ar{Aa) are the retarded and advanced vacuum propagators which for bosons 
and fermions are given by 

1 p fo + m 



(^) = i^ ^^-T' ^RiA)ik)= \+ ., ■ (n.3) 
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In the expressions for A]^(^ji^ (k) the usual limit e — > is considered, whenever this limit is defined, e.g. 

Af (A:) - Af (fc) = -2TTie{ko)5{k^ - m^), (II.4) 

where e(fco) = 0{ko) — 6'(— /cq) is the sign function. 

To account for noncquilibrium effects we follow the approxi matio ns described in detail in EG], and more 



recently in ||22| . It amounts only to replace the thermal n(ko), (11.2), by their noncquilibrium counter parts, 
i.e. in general by Wigner distributions n(kQ,X). This corresponds in terms of the cumulant expansion to 
approximate noncquilibrium correlations by the second cumulant only. As a further approximation we suppress 
the possible dependence on the center-of-mass coordinate X, essentially assuming a homogenous and isotropic 
medium. 

Practically as noncquilibrium distributions for bosons and fermions we take, 

n{ko,X)^.> n{ko,X)=^ for (II.5 

\^l-nF[\ko\)) \^ -{1 + nB(\ko\)) kg < 

with the Jiittner parameterizations 

^BiF){\ko\) ^ ,(|feo|4/T^l = el^ol/T^A ' ^"-^^ 

or for the actual estimates just the modified equilibrium distributions 

nB(F)(|fco|)= (II.7) 

expressed by introducing the fugacity parameter p3| A = e'^/"^, which is assumed to be energy independent. 
Obviously A 7^ 1 in case of (chemical) noncquilibrium. 

We have thus seen that if the separation of the macroscopic scale X from a fast microscopic one in the sense of 
the Wigner transform is valid, then the relevant noncquilibrium bare propagator has formally the same structure 
as the corresponding one in equilibrium. This, however, is not in general true for dressed propagators where 
the appearance of the additional out-of-equilibrium pinch singular term is to be expected . 

To show the origin and the possible appearance of pinch terms, let us consider for illustration, the selfenergy 
one-loop corrections to the (12) component of the scalar propagator. Denoting the selfenergy insertion by the 
matrix Ilab{p) we obtain the improved (12)-propagator using the Dyson equation, 

a,b=l,2 

= iDi2ip)+iSDi2{p)- (11.8) 



With the propagators specified in (II.l , II.3| ) and with the relations 



Un{p) = -U*2ip), ImHii = ^(ni2+n2i), (II.9) 



also valid out-of-equilibrium, we obtain 

SD,2ip) = n{po,X){Al{p) - A\{p))ReIln{p) 

+ in(po, X)iAUp) + A\{p))[nM - n2i(p)] 

+ AR{p)AA{p)[n{po,X)Il2iip)-il + n{po,X))U,2ip)]. (ILIO) 
Here the ill-defined product A/j(p)A^(p) appears, giving rise to possible unpleasant pinch singularities |T^Jl9[]. 



In case of equilibrium, however, it is well known that the last term in (11.10) vanishes due to the detailed balance 
condition, 

n2i(p)-eP°/^ni2(p). (ii.u) 
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As a consequence (11.10) simplifies to 



2me{pn)n{pn) 



(5'(p2)Ren(p) 



Inin(p) 



(11.12) 



where 



Ren(p) = Renii(p), 



Imn = -e po) I . : 
2 n(po) 



(11.13) 



expressed in terms of the equilibrium distribution n{po), and in terms of the (retarded) selfenergy Tl{po + iepo,p} 
p3|. S' denotes the derivative of the (5-function, and P the principal value. 



For noncquilibrium distributions (II. 5), for which obviously (1 + n)/n cxp(po/T), the last term in (11.10) 
does not cancel and may give rise to pinch singularities, which, however, have to be regularized p8[ | by proper 
resummation. 

The appearance of the pinch term in the particle propagator implies a non-trivial modification of different 
physical processes which are relevant in out-of-equilibrium systems. In the following we discuss how the nonequi- 
librium effects modify dilepton and photon rates in strongly interacting matter. In particular we show that 
the pinch (singular) term can be successfully regularized leading to finite rates for both dilepton and photon 
production from a hot mesonic and quark-gluon medium, respectively. 



III. DILEPTONS FROM A NONEQUILIBRIUM MESONIC MEDIUM 

As the first example of the application of the real time formalism to describe noncquilibrium effects we 
consider in this section dilepton production from a mesonic medium. 
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FIG. 1. Particle multiplicity per fm for a baryon chemical potential ^ib ~ O.lSGeV as a function of temperature T 



The absence of chemical equilibrium in a hadronic medium can be effectively taken into account modifying 
the particle distribution function fl^ by the fugacity parameter A as explained in the last section. In addition 
we simplify the distribution adopting the Boltzmann approximation nik^^X^ ~ Ae"'*^"'/-^, and assume that 
(5A = A — 1 is small (5A ^ 1 , i.e. a situation not far out-off-equilibrium. 

Thermal dileptons are mostly produced from hadronic decays, particle bremsstrahlung and hadron-hadron 
scattering. In Fig. 1 we show different particle number densities as a function of temperature for a fixed value 
of the baryonic chemical potential /i^ = O.lSGeV. As one can see in Fig. 1 pious and rho mesons are the most 
abundantly present particles in a medium in the temperature range 0.14 < T < 0.20GeV relevant for heavy ion 
collisions at SPS energy. Thus, discussing soft dilepton production in a noncquilibrium hadronic medium we 
consider first dileptons originating from tt — p scattering. 

The interaction of the charged pions with the neutral massive rho meson field and the electromagnetic 
potential is described by the Lagrangian |^ 
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(IILl) 



where = d^ — ieA^ — igpP^i is the covariant derivative, <I> is the complex pion field in the following assumed 
to be massless, is the rho and i^^jy is the photon field strength tensor. As is well known, the /Otttt coupling 
is rather large gp/'iir ~ 2.9; nevertheless we attempt an eiTective perturbative treatment. 




FIG. 2. One-loop diagram 



The thermal emission rate of heavy photons with invariant mass M, energy E and momentum q can be 
obtained from the trace of the photon selfenergy tensor H'^^ as follows [p|,p7[ 

" r*n?2(go,g). (III.2) 



dAPd^q/E 487r4M2 



The virtual photon selfenergy is usually approximated by carrying out a loop expansion to some finite order. 

) and using the closed-time-path formalism |p^-^ we calculate the dilepton rate in 
a thermal pionic medium at the two-loop level. Typical diagrams are shown in Figs. 2,3 and 4. The two tadpole 
diagra ms of 0{gj,) are not included, since they do not contribute to the discontinuity of the photon selfenergy 

ill. 



From the Lagrangian ([II.l 




FIG. 3. Two-loop diagram: selfenergy insertion. The labels a,b =1,2 denote the type ofirp vertex. The solid line with 
momentum label k corresponds to the p 



On the one loop level, Fig. 2, the photon selfenergy is expressed by the tree level particle propagators. 
Thus, the only difference with equilibrium calculations enters through the distribution function. In a rather 
straightforward way one recovers the standard expression for the Born term, tt+tt^ ^* ^ e+e~, 1^, which 
under the Boltzmann approximation and for M < nip reads: 

'^^ (1 + 2<5A)-^|F.(M)p eM-E/T) (III.3) 



dM^^ 967r 



where is the electromagnetic pion form factor and E'^ = q^ + AP. 

The only difference of the above result with the one derived in the equilibrium medium is the appearance 
of the (1 + 2(5A) factor which accounts for nonc quilibrium effects. The more complete expression for the Born 
contribution with quantum statistics (c.f. (|ll.5| )) and for arbitrary heavy photon momentum is, however, also 
known in the literature (see [Q). 

For soft dileptons M < rup we consider in the following the mass distribution per unit space-time volume 
dN/dM^d^x. fill To go beyond the one-loop approximation and to include the contributions due to tt^ — 
scattering one needs to calculate the two loop diagrams shown in Figs. 3-4. There are two types of diagrams 
contributing to the thermal dilepton rate from n — p interactions. These are the diagrams with real and virtual 
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p° vector mesons. The processes involving real p°'s are due to tttt pj*, irp — + nj* and p ~* ttttj* reactions. 
These contributions are obtained by cutting the two-loop diagrams (Figs. 3 and 4) such that the p is put on-mass 
shall. The virtual rho contributions, obtained by cutting the diagrams shown in Figs. 3-4 lead to 0{g^) order 
corrections to the Born term. 




FIG. 4. Two-loop diagrams: vertex type corrections 
In the following we will discuss real and virtual rho contribution separately. 



A. Dilepton rate with on-shell rho contribution 



We start to discuss explicitly the contribution from the diagram shown in Fig. 3, where the selfenergy 
correction dDi2{p) is inserted into the pionline with momentum p. The expression for 6Di2{p) is derived in 
Section 2 using the Dyson equation. 

Due to the structure of the (one- loop) selfenergy correctio n to t he pion propagator ( ILlOj ) a more careful 
derivation is required. Using the second and the third term of (11.10) we calculate the "real" correction to 1172, 



d'^p 



[p+p'f 



v{-fn{po,x) (zni2(p)-*n2i(p)) 



+ 



1 



(n(po,^)«n2i(p) - {l + n{po,X))mM)\^D2i{p), 



(IIL4) 



where a correspon ding t erm with [p < 
to TTp interactions (III.l) is given by 



p') has to be added. The pion selfenergy iTii2{p) at one-loop order due 



«ni2(p) = 5p 



(2^ 



{2p + kY[-g, 



{2p + kY Di2{p + k)D2i{k). 



(IIL5) 



In the second term in (III. 4) one can recognize the contribution of the pinch term proportional to the product 
An A A, which is not cancelled in case of nonequilibrium distributions ||l7|| . To regularize [Q the pinch singularity 
we have explicitly included a finite damping rate 7 of the pion. 

The pion damping rate 7 is determined by the "pole" (in the lower energy half-plane) of the retarded prop- 
agator. In the following we neglect thermal corrections due to the real part of the (retarded) pion selfenergy, 
such that 



Imn(po,p) 



Po 



(III.6) 
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evaluated on-shellp^ — for positive pion energy po- In the one-loop approximation the dominant contribution 
to the absorptive part of the pion selfenergy comes from the pions thermal distribution. In the Boltzmann 
approximation it leads to 



Imn(po = P,p) 



9 9 



n{E^,X)dE^ 



-A 



5^ m^T 

-J- t. g 4pT 

'An ip 



(111.7) 



Imll vanishes for p = 0, and it has its maximum near p ~ j^. We note that (111.7) gives a positive damping 
rate 7, as required. In order not to overestimate the contributions arising from the A^ Ayi terms we take for 
the numerical estimates a momentum independent value, namely the one at the maximum, 



.9^ 1 

Po7 - -Inin(po ^P,p) ^ A-^-T^. 

47r e 



{IIU 



We discuss further the estimate for the dominant contribution which is due to np nj* in the limit 7 — s- 0. 
In the Boltzmann approximation the contribution arising from (III. 4) to the dilepton spectrum reads [filj 



dNP'' 



2m? +M 



7-0 _ «n5|/4^) /ttTS 




du I 




du f 



dt 



(III.9) 



When using in (III. 9) the momentum independent value for the damping 7 ( Ill.q ), regularizing pinch singularities 
at 



^ — 0, the integrations can be explicitly performed. The leading term with 7 

5\ a\gl/A^) ' 



indicating the pinch singularity, 



dM^d^a 



A 



247r3 



2ml 5?/(47re)r 



ariz 

is proportional to I/7, 
(III.IO) 



which is actually independent of the coupling g^. 

To obtain the overall ra te for dilepton production with real (neutral) rho contributions one still needs to 
calculate the first term in ( 111.4 ). According to the calculations presented in in the limit of soft dilepton 
mass M < nip, this term reads 




ds e 



-■nip IT 



du 



(III.ll) 



where 



= i (m^ + m2 - s) (-) i y^(s-(mp + M)2)(s-(mp-A/)2), 



(-) 



mlMyt, Umin = ml + - t, 



(III.12) 



with s + t + u ■ 



+ M^. The u-integration may still be performed analytically. 

In the above formula we have included the corresponding contrib utions to the real processes from the 
two- loop vertex type diagrams of Fig. 4. The first integral in ( 111.12 ) corresponds to tttt — + p7*, the second to 
p —^ TTTTj* and the third to np — > ttj* process. 

Fin ally th e total nonequilibrium rate from real rho processes is obtained as a sum of ( III. 12 ) and the pinch 
term ( [III.10| ). 
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B. Higher order virtual corrections 



For heavy photon production at order 0{gp) we have, in order to correct the Born rate, to include virtual p 
contributions, which arise from the processes shown in Figs. 3 and 4 by cutting the diagrams in the proper way 
only through pion lines, without cutting the p line. 

In so me det ail we describe the contribution due to the selfenergy diagram (Fig. 3). Here we take the first 
term in (ILK) proportional to the real part of the pion selfenergy Il{p) and evaluate 



(2^ 



{p+p'f e{po) n{po, X) S'{p^) Ren(p) iD2i{p') 



(III.13) 



The temperature-independent part of Rell is absorbed in the definition of the T = pion mass (which we 
take approximately as TOtt = 0). The temperature dependent part in the one- loop case under consideration 
is weighted by the thermal distribution either for the p meson or for the pion (see |^^). The first case is of 
0(exp(— mp/T)), i.e. negligible, the second one therefore dominates and is expected to be of 0(T^/mp), due 
to the presence of the T — p— propagator in the loop (Fig. 3). 

In the following we present the result for the virtual rho contribution obtained from t he selfenergy diagram 
Fig. 3. In the limit nip ^ T, having in mind the mass region M < nip we get for (111.13) 



i6n 



'y,virtSE 
12 



(M,g = 0) 



OTT'^ 47r m,. 



(III. 14) 



In an analogous treatment the virtual T~ dependent contributions from the vertex type diagrams, namely 
from the first two of Fig. 4 have been calculated in . The final results for pion annihilation including virtual 
rho contributions may be estimated as. 



dN 



Born+virtual 



dN 



Born 



dM^d^a 



dAPdH 



7 g; 



T 



1-A- ^ ( 
TT 47r m 



(III.15) 



valid for nip ^ T, and for M < nip. It is interesting to note that the 0(g^p) T-dependent corrections ( 111.15 ) 
are negative and could be large. It suggests to perform resummations. 



C. Dileptons from a nonequilibrium medium - quantitative discussion 



In the last sections we have shown that the real time formalism of finite temperature field theory is an 
adequate tool to consistently describe dilepton production from n — p scattering in an off-equilibrium mesonic 
medium. In the following we discuss quantitative properties of the dilepton rate with the emphasis on the 
importance of nonequilibrium effects. 




FIG. 5. Dilepton rate as a function of M at fixed r=0.15GeV. The solid line represents the sum of all contributions 
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In Fig. 5 we calculate the over all dilepton rate from ir—p scattering in an equilibrium medium that is for A = 1. 
We indicate in Fig. 5 the contributions from different processes originating from the two-loop approximation of 
the virtual photon selfenergy. As expected the it~^'k~ p'^e^e~ reaction is a dominating source of very soft 
dielectrons with invariant mass M < O.SGeV. Thus, it could be a competing process with pion bremsstrahlung 
and Dalitz decay. At higher invariant mass M > O.SGeV the reactions with the on-shell rho meson in the initial 
state vrp —> TT e+e^ and p ttttj* are dominating processes calculated at the two-loop level. However, these 
contributions are still much smaller than the Born term due to pion annihilation tt+tt" — > p ^ 7* obtained with 
the effective vertex from VDM. This is mostly because of the pion form factor which near the rho resonance 
increases the rate by more than an order of magnitude. The destructive interference of the amplitudes in the 
different channels and also smaller phase space are the additional reasons for this behaviour. 

The contribution of the virtual rho meson calculated at e^g^ order should be considered as finite-T corrections 



to the Born tt+tt 
rather large. 
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e~^e annihilation process. As seen in Fig. 5 these corrections are negative and 



In order to analyse the nonequilibrium effects on dilepton production from the mesonic medium we plot in 
Fig. 6 the total rate characterized by different values of the chemical potential p. We observe that for positive 
p the rate is increasing with increasing p. However, naively one would expect an increase by a factor of A^, 
whereas the results in Fig. 6 show a much lower enhancement. Changing p from p ^ to p = 100 MeV only an 
effective increase by a factor 2 results in Fig. 6, contrary to a factor 4 exp ected f rom A^. This is mostly because 
of the negative contribution of the pinch — singular term summarized in (111.10). For negative values of p, that 
is for A < 1 the pinch term (III.IC) is positive, thus it increases the contribution of the irp — > ttj* process in the 
soft part of the dilepton spectrum. Consequently for A < 1 the spectrum is much flatter, with a less pronounced 
rho meson peak, as seen in Fig. 6. From this we deduce the importance of taking into account the non-trivial 
pinch term ( III.1C| ), which is due to the structure of the pion propagator in a nonequilibrium medium. 

The possible phenomenological implications of the above results, in particular the role of tt — p scattering 
in the description of the recent CERN experimental data on soft dilepton production, have been discussed in 

iTHJir 
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FIG. 6. Dilepton rate due to vr — p interactions for different values of the nonequilibrium chemical potential 



IV. HARD THERMAL PHOTONS FROM NONEQUILIBRIUM QUARK-GLUON PLASMA 



In this section we present the second example of applications of the real time formalism to describe the 
production of hard photons from an off-equilibrium quark-gluon plasma (QGP). 

The photon yield is one of the possible observables to signal a deconfined medium produced in ultrarelativistic 
heavy ion collisions. To lowest order in strong interactions hard photons in a QGP are predominantly produced 
from annihilation and Compton processes, 

q + q^ g + 1, q{q) + g ^ q{q) + i- (iv.i) 
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p' 

FIG. 7. The photon polarization tensor ni2 to first order in Us for the real photon production. The momentum 
labels define the notation used in the text. External vertices are to be of type 1 and 2, respectively, internal vertices are 
to be summed over a,h — 1,2 



It is thus clear that the photon rate is very sensitive to the momentum distribution of quarks and gluons. 
Recently different QCD motivated phenomenological models 16 1 suggest that the partonic medium created 



initially in heavy ion collisions reach local thermal equilibrium after very short time, ho weve r, being far away 
from chemical saturation. Thus, similarly as we have done for the mesonic medium (see (11.5), we assume that 
parton distributions can be approximated by thermal phase space distributions, however, with nonequilibrium 
fugacities. 



The rate of photon emission can be derived (c.f. III. 2) 



E 



dR 



<Pq 2(27r)' 



rni2i:(g), 



(1V.2) 



from the trace of the (12)-element Iii2 of the photon-polarization tensor calculated to leading order in the strong 
interaction of quarks and gluons. The relevant diagram in fixed order perturbation theory is shown in Fig. 7. 
The relevant (12) and (21) propagator components depend on the (off-equilibrium) distributions of quarks and 
gluons. 



For real photons the selfenergy Yii2 in (IV. 2) receives contributions only from spacelike loop-momenta < 0. 
The fixed order result from the diagram Fig. 7 turns out to contain IR singular contributions from small — > 0. 
For the equilibrium case it has been shown that this unphysical behavior can be cured taking hard thermal loop 
(HTL) contributions |2^ from all orders of the quark selfenergy consistently into account |2^,0 . In the spirit 
of the HTL-resummation program it is important to distinguish regions of hard from those of soft momenta 
[|3l[|. We choose to separate these scales along the line = —k'^ with kf. to be chosen on an intermediate scale: 



' Phard ' 



> kl ^ gT^ > -pL/t 



(IV.3) 



For soft momenta resummation of leading HTL-contributions from one-loop selfenergy insertions is accomplished 
by substituting the effective resummed propagator as indicated by a blob in Fig. 8. In an equilibrium medium 
the resulting partial rate is IR finite and is shown to exactly match the hard contribution thereby introducing 
the thermal mass parameter ruq as an effective IR cut-off. 

Turning to the off-equilibrium situation now our first task is to determine the appropriate generalization of 
the effective quark propagator on the hard and the soft scale. 



A. Quark propagator - hard momentum scale 



For hard momentum —p^ > k^ only the one-loop selfenergy correction to the quark propagation is to be 
considered. Denotin g the quark selfenergy insertion by S the quark propagator can be obtained from the Dyson 
equation similar to (II.S), but in terms of the quark distribution ng(po, X), 
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dSi2ip) = -n,ipo,X){Al- A\) i>Rell{p) i> ^ -n,{p^,X){Al + A\) ^(Si2 - S21) i> 

-ArAa - n,(po, ^))Si2 + n,(po, X)^2i] i>- (IV.4) 



The first term in this expression corresponds to virtual corrections and vanishes for real photon production. 
The third genuine off-equilibrium contribution, give rise to pinch singularity. However, due to the restricted 
kinematics —p^ > k"^ the singularity at = is never crossed so that in this region the pinch combination is 
equivalent to the principal value appearing in the second term: 

i (A^j + Ai) = P (^^^ ^ ArAa for -p2>fc2>o. (IV.5) 



With this relation all relevant terms in Eq. (IV.4) proportional to the fermion distribution n(po, X) are seen to 



cancel leaving behind the one loop corrected propagator 

1 



(IV.6) 



It has the same form as in the equilibrium case with the only difference that the quark selfenergy should be 
calculated with off-equilibrium distribution function. 



B. Quark propagator - soft momentum scale 



On the soft momentum scale, contrary to the discussion in the last subsection, we can not restrict the 
calculation of the quark propagator to the one-loop level. Here as in the equilibrium medium, one needs to 
carry out resummations to derive the effective propagator. In the nonequilibrium case, however, the analogous 
resummation program is complicated by the appearance of the extra terms as we have already seen. 




FIG. 8. The photon polarization tensor for soft momentum p ~ gT. The blob indicates the HTL-resummed 
quark propagator 



The resummed propagator for the off-equilibrium situation was discussed in ]17| , |18[ for the scalar case with 
the objective of providing a self consistent cut-off for the pinch contributions. Generalizing to the fermionic 
case resumming successive (one-loop) self-energy insertions 



may be rearranged into the effective propagator 133] , Q 



iSl2{p) = -nq{po,X){ 



-) 



' — E + iepa 



+ iepo ^ - S* - iepo 
[(1 - n,(po,^))(-«Si2) + ri,(po,^)(-»S2i)] 



- S* - iepo 



(IV.7) 



(IV.. 
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making use of the identities 

Sn = -E52, /mS = -^(Si2-I]2i), i?eE = i?eSii, (IV.9) 

where T,* denotes complex conjugation and ng(^g^{po, X) denotes the quark (gluon) nonequilibrium distribution 
functions. The effective propagator is weU defined as it is regularized by the complex retarded quark selfenergy 



The expression in square brackets of ( |IV.9| ) 



(1 - Tigipo, X))Si2 + ugipo, X)E2i = (1 - 2ng{po, + E+, (IV.IO) 

contains both leading (in the strong coupling g) 



-(E12-S21)- I^CFy EdE{ng{E,X)+ng{E,X)) = ml{\,Xg), (IV.ll) 
mliK \) = ^^Cf (x, + ^]^ Oig'T% (IV.12) 



and subleading 



1 



E+ = -(S12 + S21) --Po.9'y^ dkki—ng{k,X))il + 2ng{k,X))^Oig''T^) (IV.13) 

contributions. In the spirit of the HTL-resummation scheme keeping only leading contributions, the resummed 
propagator can be rearranged into 

where pinch-like terms have disappeared. P3,B3l 



The only difference of the result in (IV. 14) and the equilibrium resummed prop agato r resides in the one-loop 



quark selfenergy which has to be calculated with modified distribution functions (II. 5). 



C. Hard photon rate 



The la st step to derive the finite hard phot on ra te in off-e quilibr ium QGP is to calculate the photon selfenergy 
in ( IV. 2 ) with the propagators derived in ( IV. 6 ) and in ( [V.14 ). The calculations are well described in the 
literature for an equilibrium medium |30|. Here we quote only the results obtained in p^,Q. 



The partial rate from the soft part 



E. 



dR 
d?^ 



3a 



soft 



(IV. 15) 



is found to be proportional to m^, which also cuts the logarithm from below as expected. The parameter 
~ \/gT cutting the logarithm from above results from restricting the resummed calculation to soft exchanged 
momenta —p^ < k^. 



The hard part calculated with the propagator (IV.6) contains both the singular - cut-off dependent and 
regular contributions. The singular part of the rate in the off-equilibrium situation is obtained as 



E^ 



d^q 



2 3a 



A,m^(A„A,)e-^^/^ln(l||^ 



(IV.16) 
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the regular one as 



hard 



9 2aas 
e„ — 



(IV. 17) 



with 



C{E^, T, Xq, Xg) — Xq 



-i+(i-^h+(i-f;)i.s+c- 



+ AgAg 
+ Xg 



i-y + (Y-2)(7 + ln|^) + ^C'(2) + ^ln2 + (a-C-) 



2 + (i-y)7 + (i-y)ln^-C. 



(IV.18) 



The appearing symbols are Euler's constant 7, the derivative of the Riemann ^-function evaluated at C'(2) and 
the sums 



-0.04. 



(IV. 19) 



n=2 



n=2 



The result Eq.( jV.l? ) is successfully checked for consistency by noting that the constants present in the equi- 
librium result ||30|| are reproduced for = Ag = 1. 



Ca lculati ng the complete emission rate for real photons as the sum of the two partial rates Eqs.(IV.15,IV.16) 
and (IIV.ITD , 



dR 



dR 



dR"""-! 



soft 



d^q 



dR" 



hard 



d^q 



(IV.20) 



hard 



We see that the dependence on the arbitrarily chosen kc cancels and the final result for the hard photon rate 
can be written as 



X„ 



2E^T 



(IV.21) 



The above result shows that as in the equilibrium medium the generalized thermal mass is established as 
self con sistent cut-off for the logarithm ic si ngularity. The dynamical screening of the mass singularity seen in 
( IV.21 )-here given for the distributions ( II. 7 ) actually does not depend on the explicit form of the noncquilibrium 
quark and gluon distribution functions. Changing t he para meterization of the distr ibution functions enters only 
through the redefinition of the mass parameter in ([V.14) and the constant C in ( IV. IS] ), keeping at the same 
time the functional form of the rate unchanged. 



D. Photon production in an expanding quark gluon plasma 



In the following we discuss the importance of noncquilibrium effects on hard photon production. Having in 
mind the possible phenomenological implications of the results we include the space-time evolution of the off- 
equilibrium quark gluon plasma. We start with the assumption that the partonic system produced in heavy ion 
collisions achieves local thermal equilibrium after the time t^. Deviations from chemical equilibrium of quarks 
and gluons are measured initially by the fugacities A^ and A^. Beyond this point the system expands according 
to the Bjorken model [|4j and the chemical equilibration, that is the time dependence of fugacities is governed 
by rate equations. 

If one assumes the off-equilibrium plasma to undergo only longitudinal expansion then according to the 
hydrodynamical equation the time r-dependance of temperature is governed by the well known Bjorken's relation 
|34| ] , (valid for a free gas) 
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er^l'^ ^ const, (IV.22) 
where for the energy density of a nonequiUbrium plasma we take 



e=C-^X, + ^^X,)T\ (IV.23) 



where Uf 2.5 is the number of dynamical quark flavours. 

For the approach to chemical equilibrium we include only the dominant reaction mechanism for parton 
production. These are just four processes 

gg<-> ggg , gg < - > gg. (IV.24) 

With the above processes the rate equation can be written as the following master equations in terms of the 
weighted reaction rates i?2,3 fS^: 



1 dXq 

\q dr 



+ 



3 dT 
f~d^ 



9^/ 



describing together with (IV.22) the time evolution of the parameters T(r), Ag(T) and Ag(T). 




LHC 

RHIC 

RHIC equilibrium 



10 - [GeVj 

I I I I ^ 

DO 10 20 30 40 50 

FIG. 9. Distribution of thermal photons for j/ = from off-equilibrium QGP for the conditions at RHIC and LHC 
energies. The equilibrium result for RHIC is shown by the dotted line 



In Fig. 9 we calculate the space-time integrated nonequilibrium photon rate ( IV. 21 ). The initial conditions 
have been fixed for RHIC: n ^ 0.25fm, T, - 0.67GeV, - 0.064, ~ 0.34 and LHC: n 0.25fm, T, ~ 
imGeV, A* - 0.082, A^ - 0.43 as suggested by the Self-Screened Parton Cascade model [Q. 

To emphasize the nonequilibrium effects we also show in Fig. 9 the photon rate for the RHIC experiment 
calculated with equilibrium distributions of quarks and gluons. The initial temperature Ti ^ 0.43GeV is fixed 
requiring the initial energy density ~61.4GeV/fm'^. The results in Fig. 9 indicate the differences between 
equilibrium and off-equilibrium photon yields. As expected the higher initial temperature for the chemical 
off-equilibrium scenario increases the yield of large pt > 3GeV photons. 



Turning finally to comparisons of our results with previous studies we note first that the rate (IV.21) differs 
in details from the results of both Refs. and p^ , |38[ |; but indeed justifies the cut-off prescription employed 
therein. The discrepancy with respect to |[37[ can be attributed to the fact that distributions for the incoming 
particles have been treated in Boltzmann approximation. We find, however, that it is necessary to keep track 
of full quantum statistics in order for the screening to be complete. 

The photon yield calculated in Fig. 9 also differs from the results given in |^^. This is related to the new 
structure of the off-equilibrium photon rate derived in our approach. However, although our estimate shows 
a different functional dependence of the photon rate on the fugacity parameters, nevertheless the numerical 
differences with previous estimates are small being of the order of 20%. 



14 



V. SUMMARY 



We have calculated the thermal production rate for soft dielectrons produced in ofF-equilibrium pionic medium 
including relevant reactions from tt — p interactions. The calculations include the 0{gp) corrections arising from 
the two-loop contributions to the virtual photon selfenergy. We have discussed the importance of the pinch 
term which arises from the structure of the nonequilibrium effective particle propagator. The appearance of the 
pinch term modifies substantially the properties of the photon rate leading e.g. to the enhancement of the soft 
dilepton production in an undersaturated medium. 

The second important result of our analysis is that we show explicitly the absence of any additional contri- 
butions from possible pinch (singular) terms to the off-equilibrium photon production rate in a quark gluon 
plasma. As far as the soft momentum scale is concerned these terms are shown to be subleading with respect 
to the dominant hard thermal loop contributions. On the hard scale on the other hand the absence of pinch 
singular contributions is due to the restricted kinematics of real photon production. This result does not apply 
in general, e.g., for virtual photon production the relevant region of phase space is sufficiently enlarged for 
these terms to become substantial. This has been shown to hold in our estimate of dilepton production from a 
mesonic medium pT|] and it is also valid in the case of the quark gluon plasma p2]. 
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